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Abstract 

The process AB — » quarkonium +7 + X, for a heavy quarkonium state such as 
J/ip or T, is considered. The ratio of the cross-sections for this process in pp and pp 
collisions is shown to be a very sensitive probe of models of quarkonium formation. 
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Over the years, there has been considerable interest in the study of heavy quarkonia 
because the leptonic decay channels of these resonances provide the cleanest signals of 
heavy quark formation. The production of quarkonia takes place through subprocesses 
that are dominated by gluons. These processes have, therefore, been used to study gluon 
distributions in the hadron. 

The cross-section for the production of a QQ bound state factorises into two parts : 
a short-distance part, which corresponds to the production of a heavy quark pair from 
a hard collision of the incident particles, while the other part specifies how the QQ pair 
produced in the collision forms a quarkonium bound state. The short-distance part is, 
of course, computable in perturbative QCD but the formation of the bound state from 
the QQ pair can be specified only in the context of some phenomenological model of 
hadronisation. In particular, it is not known whether the hard scattering subprocess is 
sensitive to the quantum numbers of the bound state, or whether this information is 
completely screened by the subsequent hadronisation. 

For several years now, two models of quarkonium formation have been used to describe 
J/i/j and T production in high energy collisions : the colour- singlet model Q and the 
semi-local duality model 0. The latter model has also been referred to as the colour- 
evaporation model in the literature. In the colour- singlet model, one starts with the full 
QQ production amplitude and then projects out the state with the proper spin, parity 
and charge-conjugation assignments to describe a 2S+1 Lj quarkonium state. This state 
is required to be a colour-singlet, which is achieved by the radiation of a hard parton in 
the final state. The matrix element so obtained is then convoluted with the modulus- 
squared of the wave-function at the origin, |-Ro| 2 , (or its derivative, in the case of P- 
wave quarkonia), which finally yields the matrix element for the quarkonium production 
process. The wave-function and its derivative are fixed from the leptonic or hadronic decay 
widths of the quarkonia. In this model, because the amplitude for the production of a 
state with definite J PC are computed, it is possible to predict the production cross section 
for different resonances in a given family of quarkonia. It is also important to note that the 
hard scattering vertex is governed by selection rules that involve the quantum numbers of 
the quarkonium. For example in leptoproduction of charmonia, the 3 Si state (the J/ip) is 
produced in the subprocess — > 3 Sig, but the production of the corresponding P-states 
(the x's) is disallowed. 

In the semi-local duality model, the hard scattering vertex is completely blind to the 
quantum numbers of the quarkonium resonance. The quantum numbers of the QQ pair 
produced in the hard scattering are left unspecified, and the QQ pair is not required to 
be a colour singlet. The colour non-singlet QQ state is assumed to form the physical 
colour-singlet quarkonium state by multiple soft-gluon emission. In practice, one simply 
computes the open QQ production cross-section and integrates this cross-section between 
threshold (= 2mg) and the open QQ production threshold. Arguments based on semi- 
local duality are then invoked to relate this integral of the open QQ production cross- 
section to the sum of the resonance production cross-sections. It is clear that the cross- 
section for producing a particular resonance with a given J PC cannot be predicted in this 
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model. The cross-section for the production of a specific bound state is a fraction / of the 
integral of the open-charm production cross-section over the mass range considered. This 
fraction / is, however, not a prediction of the model but can be obtained by comparison 
with experiment. 

The colour-singlet model has been applied to the study of quarkonium production in 
both lepton-hadron and hadron-hadron collisions 0. In leptoproduction, the heavy- 
quark pair is produced through the fusion of the incident virtual photon with a gluon 
from the nucleon. The predictions of this model have been compared with the data on 
J /if) production at large p T and small z in lepton- nucleon collisions 0, |J, where px is 
the transverse momentum of the J /if) and z = Ej/^/E^; Ejj^ and E y are the energies 
of the J/ip and the incident virtual photon, respectively. The model gives an adequate 
description of the pt and rapidity distributions of the J /if) for large px and for z < 0.7 - 
0.8. The overall normalisation predicted by the model turns out to be lower than the 
data by a factor of 2.3 M. Similarly, data on J/ if) and T hadropro duct ion at large px also 
seem to be consistent with the predictions of this model |J, except for the normalisation. 
Comparisons of the predictions of the semi-local duality model have also been made with 
data on leptoproduction and hadroproduction of J /if) 0. The remarkable thing is that 
this simplistic model is also capable of explaining the data at a level comparable to 
the more elaborate colour-singlet model. A recent analysis of J /if) production in // - p 
scattering suggests that the colour-singlet model is somewhat favoured in the region 
z < 0.8 as compared to the semi-local duality model, but the experimental errors are too 
large for any definite conclusions to be drawn. 

In this letter, we suggest that the production of a photon in association with the 
quarkonium will help to clearly distinguish between these two models of quarkonium 
formation. The associated production of J / if) and photon at large px has been studied in 
the context of collider ||, fixed-target |J and polarised pp scattering experiments [fT0[] , in 
the framework of the colour-singlet model; it has been shown to be a very sensitive probe 
of gluon distributions. This is so because, in the colour-singlet model, the only subprocess 
that contributes to J /if) + 7 production, with the J / if) and 7 produced back-to-back at 
large px, is the following: 

g + g^J/iP + 1 . (1) 
The colour-singlet model predicts that no P-states are produced because the subprocess 

9 + 9 - X + 7 (2) 
is disallowed because of C-invariance. Moreover, the subprocess 

9 + 9 ~^X^ J/ip + r T (3) 

does not produce a photon at large px and does not contribute to J /if) + 7 production. 
Finally, the gg-initiated subprocess 

q + q -> J /if) + 7 (4) 
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does not contribute, because the cc thus produced will not be in a colour-singlet state. 

On the other hand, the gg-initiated subprocess will contribute in the semi-local du- 
ality model. This indicates that these models will strongly differ in their predictions for 
quarkonium+7 production. This difference can be particularly enhanced if we study the 
ratio of the cross-sections for quarkonium+7 production in pp and pp collisions. In what 
follows, we study the ratio 

daPP/dp T 

In the colour-singlet model, the ratio R is clearly equal to unity, because only the gluon- 
initiated subprocess contributes. Moreover, this will be true for all pt and will also 
be independent of the quarkonium resonance considered. Of course, this is all true at 
the lowest order in QCD : higher-order diagrams in the colour-singlet model will have 
contributions from both the gluon- and quark-initiated processes. But the ratio R, which 
we propose to study, will be more resistant to the effect of higher-order corrections than 
the absolute cross-sections themselves. So we expect that R will not deviate significantly 
from unity in the colour singlet-model, at least at fixed target energies. However, even 
at the lowest order in perturbation theory, the semi-local duality model allows for both 
gluon- and quark-initiated processes. Hence, we expect that R will differ from unity in 
this model. 

We present, in the following, the pt dependence of the ratio R in the semi-local 
duality model for J/ip + 7 production. In this model, we need to integrate the open 
charm production cross- sect ion, given by the processes 

9 + 9 -> c + c + 7, 

g + g^c + c + 7, (6) 

with the invariant mass of the cc pair restricted to lie between Am 2 c and 4m|,, the lat- 
ter being the threshold for open charm production. The kinematics is identical to that 



described in Ref. [11 1. The four- momenta of the initial partons are taken to be 



fs 

— (1, sin# cos0, sin^ sin0, cos#), 
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~2~(1j — sin6* cos0, —sin6 sin0, — cos6>). (7) 

The xy plane contains the final particles, and the direction of the photon is taken to define 
the +ve x-axis. The four-momenta of the c, c and the photon are 

fs 

Pc =— 2:3(1, As cos6» 35 , (3 3 sm6> 35 , 0), 
Pc =—£4(1, At cos6>45, At sin0 45 , 0), 



Pl =^U(1, 1, 0, 0), (8) 
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A = yJl-4m*Ms. (9) 

In the above equations, s is the subprocess centre-of-mass energy and is related to the 
hadronic centre-of-mass energy S by the usual relation, s = X1X2S, where x\ and X2 are 
the momentum fractions of the incident hadrons carried by the partons. From energy 
conservation we obtain, 

x 3 + x 4 + x 5 = 2 ( Xi < 1, i = 3,4,5). (10) 
The angles 6*35 and #45 are given by 

(to) 2 - (to) 2 - x\ 



COS0 



2f3 3 x 3 x 5 



a (to) 2 - (to) 2 - x\ 

COS ^ 45 = 2to^ • 

The pt of the photon (which is the same as that of the J /if) is given by 

y/s. 



p T = 2 X5 v cos2 ^ + sin 0sin 0. (12) 



The differential hadronic cross-section is then given by 
da aa 2 s pT 



dx\dx2dx2,dx 1 id(j)dpT A.nsx 2 cos6cos 2 



H(x 1 ,x 2 )\M qq \ 2 + G(x 1 ,x 2 )\M gg \ 2 , (13) 



where H(xi,x 2 ) and G(xi,X2) are parton distribution factors for the qq and the gg sub- 
processes, respectively. |M g q| 2 and |M 33 | 2 are the squared matrix elements for the sub- 
processes given in (||). These matrix elements are obtained from those for heavy quark 
leptoproduction by crossing. The matrix elements for the latter process are given in 
Ref. 0. 

Integrating over all variables other than p T in Eq. QI3]), we obtained the cross-sections 
for pp and pp collisions as a function of pr- The integration was done by restricting the 
invariant mass of the cc pair as described above, and also restricting the rapidities of 
the J /if) and the photon to be both equal to zero. The computations are done for fixed 
target energies at Fermilab (y/s ~ 38.75GeV). For the parton densities Owens' Set 1.1 
distributions |H| were used, and the Q 2 for the process was taken to be Am 2 + p\. 

The ratio R of the cross-sections computed in this model (curves labelled II) for J / if) 
production is shown in Fig. 1, as a function of p?- The ratio is significantly smaller than 
unity and also shows a strong p T -dependence. With increasing p T , larger values of X\ and 
X2 are sampled and the qq process becomes more important. Consequently, there is a 
larger suppression of the ratio R. We have also computed R for the case of T production, 
and find that it also depends on the mass of the quarkonium. In the case of T, we find that 
R is smaller than that obtained for J /if). In Fig. 1, we have also shown the curve (labelled 
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I) for R expected in the colour-singlet model. The ratio is simply unity, independent of 
Pt and the mass of the quarkonium. 

The experimental feasibility of measuring this process at fixed-target energies was 
dealt with in detail in Ref. 0. We would like to simply point out here that one could 
get a enhancement by a factor of 100 if a nuclear target were used and the cross-sections 
from pA and pA were compared. Nuclear effects would cancel in the ratio. 

In conclusion, we find that the production of a photon in association with a quarkonium 
such as J/ip or T, will provide important insights into the nature of quarkonium formation. 
In particular, by studying the ratio of the pp and pp induced cross-sections at fixed-target 
energies it will be possible to discriminate between the colour- singlet and the semi-local 
duality models of quarkonium formation. Just as this work was being completed, we 
recieved a preprint ]nj where leptoproduction of J/ip + 7 is discussed as a probe of 
models of quarkonium formation. 
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Figure caption 

Fig. 1 The ratio R as a function of px- I and II denote colour-singlet and semi-local duality 
models, respectively. The quarkonium states to which the curves correspond are 
shown in parenthesis. 
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